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O ; Abstract 

The proportionality between single charge-exchange reaction cross sections in the forward direc- 
ts ' tion as found, for example from (p,n) and (^He,t) and from {n,p) and ((i,^He) reactions, and the 
^ ■ Gamow- Teller (GT) strength into the same final nuclear states has been studied and/or assumed 
often in the past. Using the most physically justified theory we have at our disposal and for the 
specific example of the ''^Ge-''^Se system that may undergo double beta-decay, we demonstrate that 
I the proportionality is a relative good assumption for reactions changing a neutron into a proton, 
\^ ' i.e. ''^Ge(p, n) ''^As. In this channel, the main contribution to the GT strengths comes from the 

■ removal of a neutron from an occupied single-particle (SP) state and putting a proton into an unoc- 

■ cupied SP state having either the same state quantum numbers or those of the spin-orbit partner. 
. In contrast to this, in the second leg of the double beta decay a single proton must be taken from 

^ I an occupied SP state and a neutron placed in an unoccupied one. This second process often is 

I Pauli forbidden in medium-heavy nuclei and only can be effected if the Fermi surface is smeared 

• out. Such is the case for ^^Se(n,p)'^^As. Our results suggest that one may not always assume a 

■ proportionality between the forward-angle cross sections of the charge-exchange reactions and the 
^ . GT strength in any such medium-heavy nuclei. The discrepancy originates from a pronounced 
^ I effect of the radial dependence of the nucleon-nucleon (NN) interaction in connection with the 

I Pauli principle on the cross sections in the (n,p) reaction channel. Such a radial dependence is 

^ ' completely absent in the GT transition operator. 

I. INTRODUCTION trinoless double beta-decay transition prob- 

ability is important since it, in conjunction 

with measured data and assuming that the 
In recent years, interest in the relation be- ^.^^^ ^^^^^.^^ exchange is the leading con- 
tween charge-exchange reactions in the for- ^.ibution, defines an absolute scale for the 
ward direction and the GT strength has in- ^^^^ Majorana neutrino (for reviews 
creased due to the connection of the GT ^^^^ ^ 0,1,1). The two-neutrino double 
strength with the two-neutrino doub e beta ^^^^ ^^^^^ ^^^^.^ ^^^^^^^ ^2. g^^^^ a 
decay [i, 1^. The latter process may help to ^^^^ intermediate 1+ states of prod- 
test the nuclear wave functions required in ^ct of the GT transition matrix element from 
calculations of the matrix elements for the ^^^^^ ^^^^^ ^^^^^ .^^.^^ ^^^1^^^ 
neutrinoless double beta decay. The neu- ^^^^^^^ reQ^) ^ an 1+ state in the interme- 

diate nucleus (''^As) and the GT transition 
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divided by the corresponding energy denom- 
inator. However, a test of the two-neutrino 
double beta-decay calculations using the GT 
strengths extracted from the measured elec- 
tron capture (EC) and the single-beta decay 
of the intermediate nucleus is only possible if 
the ground state in the intermediate nucleus 
is a 1+ state (not the case for '^^As), and if the 
two-neutrino double beta decay is dominated 
by the transition through this state. Our in- 
terest with the mass-76 systems in part stems 
from the fact that for the neutrinoless dou- 
ble beta-decay matrix element in past evalu- 
ations gave values ranging from 2.23 to 5 ac- 
cording to the chosen model of nuclear struc- 
ture jfi]. 

Complementary to the direct measure- 
ment of the GT strength by the EC cap- 
ture in the first leg followed by a /3~ tran- 
sition from the lowest 1"*" state of the in- 
termediate nucleus, are single-charge trans- 
fer reactions like [p, n) and (^He,t) on the 
ground state of the initial nucleus and {n,p) 
and {d, ^He) to the ground state of the fi- 
nal nucleus; all connecting by the interme- 
diate l"*" states. If the forward cross section 
of these charge-exchange reactions are pro- 
portional to the corresponding GT strength, 
the two-neutrino double beta-decay probabil- 
ity calculations can be checked, although the 
information about the relative phases of dif- 
ferent contribution can not be extracted from 
the experimental B{GT). Then it is possible 
to test the quality of the calculations for the 
neutrinoless double beta decay. 

The proportionality between the forward 
single charge-exchange cross section and the 
GT transition probabilities has been studied 
extensively in the past. Refs. 0] and [sl are 
particular contributions. Taddeucci et al. 0] 
present an interesting analytic study of the 
proportionality involving the single charge- 
exchange reaction (p, n) cross section at zero- 
momentum transfer which we take as quite 
typical of all studies of the problem. They 
assumed that only angular momentum trans- 
fer L = is important at forward scatter- 



ing angles and that the eikonal approxima- 
tion is valid to describe the relative motion 
wave functions of the incoming and emergent 
nucleons. Under those assumptions, they ob- 
tained an expression for the proportionality 
between the forward charge-exchange reac- 
tion cross section and the GT transition prob- 
ability, both to l"*" states. However, in that 
study 0] there are a number of other as- 
sumptions, many of which are questionable. 
First, a single particle-hole configuration is 
assumed for the structure of the nuclear tran- 
sitions. They also assume that the reaction 
mechanism can be taken in either a plane or a 
distorted wave impulse approximation. Fur- 
thermore, the radial wave functions for the 
initial neutron and for the final proton in 
the {p, n) reaction are assumed identical and 
the effects of antisymmetrization between the 
projectile and the target nucleon is treated 
rather crudely. Limiting themselves to use 
the impulse approximation means that they 
use A^A^ amplitudes in calculations and not a 
specific finite ranged A^A^ interaction. The 
use expressions given by Franey and Love 
which were derived from the SP84 ampli- 
tudes for free A^A^ scattering. The nuclear 
medium has dramatic effects in making the 
effective interactions between projectile and 
every bound nucleon in the target quite dif- 



ferent to the free A^A^ case [10|. Similar con- 
cerns exist even with some of the limitations 
are removed, for example by using the dis- 
torted wave impulse approximation or by us- 
ing phenomenology to define relative motion 
wave functions in a distorted wave approxi- 
mation (DWA) approach. Such concerns we 
outline later in the text. 



Ejiri je], m has also made extensive study 
of the proportionality link. In his review |6|, 
the proportionality of the forward charge- 
exchange cross section for {p, n) and (^He,t) 
are shown in Figs. 10 and 15 for Fermi 
and GT transitions, respectively. The pro- 
portionality of the charge-exchange reaction 
cross section to the GT strength correspond- 
ing to a {n,p) reaction is depicted in that 
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review by using (rf, ^He) in Fig. 12, by using 
(t,^He) in Fig. 17 and by using ("^Li, "^Be) in 
Fig. 18. Ejiri found that the proportional- 
ity with the forward scattering cross sections 
from single charge-exchange reactions with 
type {p, n) to the GT strength was good for 
all nuclei to mass A = 124. However, pro- 
portionality studies for the charge-exchange 
reaction of the type {n,p) was investigated 
only for masses to A = 12. Nuclei rele- 
vant for the double beta-decay proportional- 
ity of those charge-exchange reactions to the 
GT strengths were not considered. Such are 
needed of course as they are important for the 
two-neutrino double beta decay in the second 
leg where a proton changes into a neutron. 
Often that change cannot be effected by the 
GT operator r'^a which can only change par- 
ticle types in orbits having the same quantum 
numbers or into the spin-orbit partner of that 
level. 

There have been many previous studies 
seeking nuclear matrix elements from exper- 
imental data with ref. [l2| the most recent. 
Cole et al. [l^ studied charge-exchange reac- 
tions from ^^Ni. This nucleus is medium mass 
but it has a neutron excess of only 2. The de- 
gree of Pauli blocking to differentiate between 
isospin raising and lowering transitions then 
is small. The case should be classed with 
those of most light mass studies. The cases 
we consider on the other hand have a sizeable 
neutron excess and so the Pauli blocking ef- 
fects in the (n,p) reactions are much more 
important than in the (j9, n) cases. 

We consider the proportionality question 
again but make use of the best available re- 
action codes to evaluate cross sections for 
the charge-exchange (p, n) and reac- 
tions. We consider specifically the very pop- 
ular double beta-decay transitions ^^Ge 
''^As — > ^^Se. The nuclear structure of the 
initial, the intermediate and the final states 
in these nuclei have been defined using the 
Quasiparticle Random Phase Approximation 
(QRPA) with realistic forces (Bonn CD po- 
tential (iSj) and with matrix elements calcu- 



lated by solving the Bethe-Goldstone equa- 



tion |14j . Il5l |. The results show that the 
forward charge-exchange cross section of the 
type (p, n) for the first leg of the double 
beta decay is nicely proportional to the GT 
strength but that the forward reaction cross 
section of type {n,p) shows rather large de- 
viations from this proportionality. The latter 
is due to Pauli blocking since a proton from 
an occupied level must be transformed into a 
neutron in an empty level with the same (n, 1) 
quantum numbers. Due to the radial de- 
pendence of the A^A^ interaction, the charge- 
exchange reaction can proceed by transition 
between SP orbits that differ in (n,£). Such 
effects violate the proportionality between 
the forward charge-exchange {{n,p) (d, ^He), 
{t, ^He), ('^Li,'^Be) . . . ) cross sections and the 
GT strength. 

Quality of the QRPA approach for de- 
scription of the GT strengths and double beta 
decay has a long history 0, i, [3, H 
liel, H [li, H [Ifl]. The nuclear wave func- 
tions calculated within the QRPA have been 
shown to provide good description of differ- 
ent properties of giant multipole resonances 
and low-lying collective 2"*" and 3~ states. 
The gross structure of the GT strength dis- 
tribution as well as the position of the GT 
resonance in the intermediate nuclei is cor- 
rectly reproduced within the QRPA provided 
that the particle-particle strength of realistic 
nucleon-nucleon interaction is slightly renor- 
malized by a factor 0.8 < gpn< 1.0, depend- 
ing on the model basis size 0,13, 0] • 

The nuclear shell model, which nicely de- 
scribes nuclear states in the sd shell of pos- 
itive parity and where it is more reliable 
than the QRPA, is not able to describe 
the states relevant for double beta decay in 
the pf and sdg shells. The Strassbourg- 
Madrid collaboration 2l|] can only handle 



a basis consisting of four single-particle lev- 
els (1/5/2,2^3/2, 2pi/2,l5'9/2) for the nuclei in 
the vicinity of ^^Ge. Since the spin-orbit 
partners I/7/2 and l(/7/2 are missing in the 
model space, the model-independent Ikeda 



3 



sum rule [22| is strongly violated and the GT 
strength calculations for ^^Ge within the shell 
model are not trustful. 

The QRPA model for the nuclear wave 
functions is considered in the next Section 
while that of charge-exchange reaction theory 
is developed in Sect. IIIII Then, in Sect. [IV] 
we present the results and give conclusions in 
Sect. M 



II. NUCLEAR WAVE FUNCTIONS 

The majority of calculations of the two- 
neutrino and the neutrinoless double beta de- 
cay have been made using the QRPA [li, 



151 . Il6l . llTj . Although the starting points 
of all these studies are very similar, ma- 
trix elements calculated for the neutrino- 
less double beta-decay transition probabili- 
ties differ. For example, those obtained in 
Refs. IJ, |l5|, |l6| are quite different from the 
ones of Ref. [l7|. This is a reason to seek 
tests of wave functions by deriving the two- 
neutrino double beta-decay probability from 
the GT strengths between the initial and fi- 
nal nucleus to a large number of 1"*" states in 
the intermediate nucleus. 

Herein we use wave functions obtained 



from QRPA calculations |14j . Il5[ | in which the 
Brueckner reaction matrix elements of the 
Bonn CD potential [13] for the A^A^ interac- 
tion were used. The strength of the NN ma- 
trix elements in the particle-particle channel 
has been slightly adjusted, by a factor gpp, 
to reproduce the experimental two-neutrino 
double beta-decay probability. For ^^Ge, this 
value is gpp = 0.85 for a 9 level basis {pf and 
sdg major shells). We used the unquenched 
values gph = 1 and gA = 1-25 for the particle- 
hole channel renormalization factor and the 
axial coupling constant gA, respectively. 



tization as 

(^JM = {pmp\bjM\nm 

pn,mpmn 

= J-'J2{p\\bj\\n) C\pn,JM).{l) 

pn 



In this equation, j = \J2J + 1, 

C'^{pn,JM) = [aj, anY^\ and bjM 
can be r~ (Fermi), ar^ (GT), or any 
other operator, including ones that have 
r-dependence. The time reversed creation 
operator is defined as a]^ = (— )-'~"^a]_^. 
Edmond's version of the Wigner-Eckart 
theorem has been used. The definition of 
the spherical harmonics includes a factor 
in order to ensure the above expression for 
the time-reversal operation. 

The reduced matrix element of such SP 
operators between the ground state of a 
mother nucleus and an excited state of the 
daughter nucleus is given by 

= J-'J2^p\\bj\\n) Q^-\pn,J),{2) 

pn 

where the elements of the transition matrices 
g^~\pn, J) are the reduced matrix elements, 

g^-\pn,J) = {r\\C\pn,J)\\0+) . (3) 

The corresponding formulae for the (3^- 
channel are obtained by the changes 

C\pn,J) ^C{pn,J) 
g^'\pn,J) g^^\pn,J). (4) 

In the RPA, a nuclear state having angular 
momentum J and projection M, is created 
by applying the phonon operator Qj^^ to the 
vacuum state |0^p^) of the initial, even-even, 
nucleus, i.e. 

\JM) = QIj,,Kp^); Q,m|OW = 0. (5) 



Any single-particle (SP) operator of the Introducing the quasiparticle creation and 
/9"-type can be represented in second quan- annihilation operators, a+ and 
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p, n) defined by the Bogolyubov transforma- 
tion, 



M ^ (6) 



tlie plionon operator Qj^// can be written 
witliin tlie QRPA as 



pn 



-Y^^U{pn,JM)], (7) 



wfiere 



A\pn,JM) = [al^ai] 



JM 



and tlie forward- and backward-going free 
variational amplitudes X and Y satisfy the 
matrix equation, 



A B 
B A 



XT. 



1 
-1 



(9) 

Here m identifies different roots of the QRPA 
equations for a given .P and 

^= (o^PAl[A[i^,^1]|oW> 

^ = 'KpA\[AdH,A]]\Qlp^). (10) 

For a realistic residual interaction, the matri- 
ces A and B are 



aJ" 



gppG{pn,p'n';J) 

X {UpUnUp'Un' + UpUnUp'Un') 

- 9ph F{pn,p'n'; J) 

X {UpVnUp'Vn' + VpUnVp'Un') 

Qpp G{pn,pn'; J) 

X {UpUnVp'Vn' + VpVnUp'Un') 

■ Qph F{pn,p'n'; J) 

X {UpVnVp'Un' + VpUnUpiVn') 



where G{pn,p'n',J) and F{pn,p'n',J) are 
particle-particle and particle-hole interaction 
matrix elements of a G-matrix, respectively. 



B, 



pn,p'n' 



Within the QRPA, one has 

C^{pn,JM) = UpVnA\pn,JM) 

+ VpUnA{pn,JM), (11) 

and the transition matrix takes the form, 

g^-\pn, J) = j{upVnXli^ + VpUnYp^;^'^), 

g'^+\pn, J) = j{vpUnXg^ + UpVnY^'l^). 

(12) 

Correspondingly, the B{GT) values for the 
GT transitions 0^ 1+ can be written as 

i?(GT(-)) = 1(1+11 5^ a„r-||0+)r 

a 

pn 

X [UpVnXp^ ^ + VpUnYp^ ^) , 

5(GTW) = 1(1+11 5^ a„r+||0+)r 

a 

= I '^HWWp) 

pn 

X {VpUnX^l^^^ + UpVnYp^^''^)\\ (13) 

In calculations, a harmonic oscillator with an 
oscillator length parameter b = 2.09 fm has 
been used to specify the SP wave functions 
for ^^Ge and ^^Se. Those functions are pos- 
itive at the origin. Using N=3 and N=4 os- 
cillator shells in the QRPA calculations for 
transitions to 1+ states in ''^As gives a set of 
23 two-quasiparticle excitations per Eq. ([8]) 
to be included, via Eqs. ([7]) and (Q, into the 
QRPA phonon creation operator for 1+ states 
in ''^As [l4, 15|. That set is shown in Tabled 
The individual components are identified by 
the label ID which will be used in the dis- 
cussion of results. Of those two-quasiparticle 
states, the ones labelled with ID = 5, 6, 15, 
16, 21 and 22 cannot be excited by the GT 
operator. 
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TABLE I: Two-quasiparticle configurations 
forming the QRPA structure of 1'*' states in ''^As 
(relative to the ground state in "^^Se) and corre- 
sponding without pairing to particle-hole states. 





q-p 


q-h 




q-p 


q-h 


ID 


n£j 


nij 


ID 


nij 


nij 


1 


0/7/2 


0/7/2 


13 


O57/2 


O59/2 


2 


0/7/2 


0/5/2 


14 


O57/2 


O57/2 


3 


0/5/2 


0/7/2 


15 


O57/2 


14/2 


4 


0/5/2 


0/5/2 


16 


14/2 


057/2 


5 


0/5/2 


lP3/2 


17 


14/2 


14/2 


6 


lP3/2 


0/5/2 


18 


14/2 


14/2 


7 


lP3/2 


lP3/2 


19 


14/2 


14/2 


8 


lP3/2 


lPl/2 


20 


14/2 


14/2 


9 


lPl/2 


lP3/2 


21 


14/2 


2Sl/2 


10 


lPl/2 


lPl/2 


22 


2Si/2 


14/2 


11 


0^9/2 


Oc/9/2 


23 


2Sl/2 


2si/2 


12 


0^9/2 


Oc/7/2 









In Fig. [H a set of one-body density ma- 
trix elements, of Eq. (OBDME here- 
after), for the excitation of five particu- 
lar 1+ states out of the total 23 found is 
shown for each of the 23 components (ID). 
These are the states of special interest regard- 
ing the '^^Se{n,p) zero-degree and/or zero- 
momentum transfer cross sections considered 
later; being the strongest of the 23 charge- 
exchange excitations considered. Of note is 
that the strongest OBDME of the fourth 
state belongs to the two-quasiparticle state 
1^3/2 — 0/5/2 (ID=6). That component, read- 
ily excited in the charge-exchange reaction, 
cannot be excited by the GT operator. Thus, 
one may anticipate that the proportionality 
between the charge-exchange reaction cross 
section and the GT strength for this state 
may be different to those of others. 

III. REACTION THEORY 

The charge-exchange reaction is described 
in a DWA in which one requires transition 




5 10 15 20 
QRPA 1 components 



FIG. 1: (Color online) One-body density ma- 
trix elements p^^\pn,l'^) of Eq. (jl2|) . in arbi- 
trary units, for six QRPA states 4, 7, 9, 10, 
12, 13 defined in Table |TI] as functions of the 
two-quasiparticle components labeled according 
to Tabled! 

structure details, optical model wave func- 
tions (the distorted waves), and a transition 
operator by which the reaction is effected. 

A. The optical potentials 

Conventionally distorted wave functions 
are the relative motion wave functions ascer- 
tained from an optical potential with which 
good fits to elastic scattering data (cross sec- 
tion and spin observables) have been ob- 
tained. In many studies those potentials have 
been assumed to be local in form and usu- 
ally of Woods-Saxon type with parameters 
adjusted to find a good fit to elastic scat- 
tering data. But the associated relative mo- 
tion wave functions are not guaranteed to be 
proper. Only the asymptotic (large radius) 
properties are tested by such data fitting 
since one only requires the scattering phase 
shifts which are then used in Legendre poly- 
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nomial sums that define the cross sections 
etc. There are other concerns about phe- 
nomenological model potentials to be noted. 
The most serious concern is the violation of 
the Pauli principle. A local potential assum- 
ing smooth parameter variations with energy 
will support bound and resonance states of 
the compound system, and the projectile can 
be captured in any. But the set includes 
states that are densely occupied in the ac- 
tual target. Further the resonances are then 
all single particle in nature and it is known 
that all nuclei support compound and quasi- 
compound ones as well. For low energy scat- 
tering one needs a better reaction theory, 
such as the multi-channel algebraic scatter- 
ing theory (MCAS) to do better. Studies 



24 have re- 



with that MCAS approach [2 
vealed that when coupled-channel effects are 
important then the effects of violation of the 
Pauli principle with potential-like models are 
most severe, and worse can lead to erroneous 
Physics. This is one reason why we consider 
energies at which specific coupled-channel ef- 
fects, such as of virtual excitation of the 
giant resonances, are minor if not negligi- 
ble. Coupling to states in the continuum 251 
has been investigated, but, without taking 



into account the effects of the Pauli princi- 
ple in the optical potentials, that result lacks 
some credibility. No doubt such may have 
a marked role in reactions such as break-up 
and inelastic scattering into the continuum, 
but we do not consider such coupling to the 
continuum as the transitions of interest are 
to specific isolated states of given spin-parity. 
For such, we have found no case as yet for 
which the basic (/-folding method (described 
next), when defined with good spectroscopy, 
requires additional reaction processes to give 
reasonable results. 

A more physical approach is to form opti- 
cal potentials by folding an A^A^ interaction 
with the target ground state structure. In 
that way one can also ensure that the Pauli 
principle is not violated. However, once the 
target structure has been set, then one has 



to choose the A^A^ interaction. For some 
time now it has been known that the interac- 
tion differs from the free A^A^ one. Medium 
effects lead to the effective A^A^ interaction 
being energy- and density-dependent as well 
as complex. The current best practice, at 
least for energies below 3-3 resonance exci- 
tation, is to use an effective interaction built 
from the A^A^ (/-matrices that are solutions of 
the Bethe-Brueckner-Goldstone (BBC) equa- 
tions \ld\. Using those (/-matrices, both on- 
and off-shell values and for 32 A^A^ angular 
momentum channels, lead to an effective A^A^ 
interaction in coordinate space that is a mix- 
ture of central, two-nucleon spin-orbit, and 
two-nucleon tensor components. Details of 
that mapping are given elsewhere l3 • 

Of great importance is that on using such 
an effective A^A^ interaction in forming op- 
tical potentials, and when account is taken 
of the Pauli principle, those optical poten- 
tials are strongly non-local and partial wave 
dependent. Non- locality arises from the al- 
lowance for knock-out (exchange) arnplitudes 
in the so-called (/-folding procedure m. Do- 
ing so requires more than just the densities 
of the nuclear ground state. One requires the 
ground state OBDME 



Pgs 



(0,+sll[al®a,](^=°)||0+). (14) 



Assuredly, the relative motion wave functions 
will differ from those found using phenomeno- 
logical (local) potentials even if the potentials 
are phase equivalent. The Perey effect is one 
ramification. Exchange effects are also most 
important in evaluations of non-elastic scat- 
tering and that will be discussed later. 

In coordinate space, the (/-folding optical 
potential can be written 

U{t, r'; E) = 6{r - r') j p(s) (/^(r, s; E) ds 
+ ^n,v^*(r)(/^^'(r,r';i?) v^,(r'); 

i 

Ms) = 5^^iV^*(s)v5i(s). (15) 
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Here p(s) is the nucleon density for nucle- 
ons with the occupancies rii. To evalu- 
ate these potentials requires specification of 
three quantities. They are the single nucleon 
bound state wave functions v^i(r), the orbit 
occupancies nj, which more properly are the 
nuclear OBDME of Eq. ([H]), and the NN 
5f-matrices g^^^^{r, s; E). 

B. The effective interaction between 
projectile and bound nucleons 

The 5f-matrices in the equation above, are 
appropriate combinations of A^A^ interactions 
in the nuclear medium for diverse A^A^ angu- 
lar momentum channels. For those A^A^ in- 
teractions, much success has been had using 
an effective A^A^ interaction, now commonly 
designated as the Melbourne force [lo|, and 
which has the form g^^ 



where r = Tq — ri and = ||(ro + ri)|. It 
is based on the ^--matrix of the Bonn B po- 
tential 261]. In the prescription, the Fermi 
momenta relate to the local density in the nu- 
cleus at distance R from the center when ri 
are the coordinates of the colliding projectile 
and bound nucleons. {ST} are the spin and 
isospin quantum numbers of the A^A^ system. 

For use in the DWBA98 program [27| . 
these effective A^A^ ^(-matrices are, specifi- 
cally. 



9'eJf = 9!Jf{r.E- kf) 



1=1 

Y.gfff{r,E-kf) e. 

i=l 



— a. 'r 

e ^ 



[S,T] 



(16) 



where Gj are the characteristic operators for 
central forces {i = 1), {1, (a ■ a), (r ■ r), (a ■ 
ar ■ r)}, for the tensor force {i = 2), {S12}, 
and for the two-body spin-orbit force (z = 3), 
{L- S}. The Sj'\E] kf) are complex, energy- 
and density-dependent strengths. The prop- 
erties of the (/-matrices are such that, not 



only can the ranges of the Yukawa form fac- 
tors be taken as independent of energy and 
density [lo|, but also four suffice with this 
approach for energies to just below the 3-3 
resonance threshold. 

The strengths (and ranges) in these effec- 
tive A^A^ interactions were found by mapping 
their double Bessel transforms to the A^A^ g- 
matrices in infinite nuclear matter (solutions 
of the BBG equations). With a: {LL'JST}, 
this mapping is 

^?effrLL'(^''^'^) = E(©^)^- (17) 

i 

where the radial integrals expand to 

J, = / r^^^jLiq'r) Q^^^^ir, E; kf) jy^qr) dr 



e 



(i) 



X / ^ jiiq'r) jL'{qr)dr 

Jo r 

J2sf\u) r"(g',g;/.f). 



(18) 



Therein A = 2 for the tensor force. In appli- 
cation, a singular valued decomposition has 
been used to effect this mapping. 

C. The DWA for non-elastic reaction 
analyses 

In the DWA, amplitudes for a non-elastic 
scattering of nucleons from nuclei, through a 
scattering angle of 9, and between the states 
|j(,Mj) and |j/,M/), are 



TnwA = rr/r (9) 

x:.7^(koO) (^^j^.Mf{l---A) 
Aj29eni^o,i,E; kf)PsPT 

ST 

Aoi{\xi'\kS))\^j.MA^---A))}, 

(19) 



X 
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where are the spin quantum number of 
the nucleon in the continuum, x^^^ ^ire the 
distorted waves, and gfj^ijQ^i, E;kf PsPt is 
the spin-isospin Melbourne force. The oper- 
ator effects the antisymmetrization of the 
two-nucleon product states. 

Then, by using cofactor expansions, 

\mjM) = A-^l'^ J2j,m l^jM ) «im|^JM), the 

matrix elements become 



rj,MfM,u'v 

JfJi 



ji,02i,S,T 



X {x[' ^(koO)|(v^,,^,(l)| (7f//ro,i, E- kf) 



(20) 

The density matrix elements in the ampli- 
tudes reduce as 



JfMf 



J2hlY''""''Hji,j2,m,,-m2\I,N) 

I{N) 

IN 



X 



JfMf 



J,Mi 



. 1) {j,,j,,m,,-m,\I,N) 

I{N) 

X {Ji, I, Mi, N\Jf, Mf) .] ^ S,,j,i, 



(21) 



where Sj^j^i are the transition OBDME. The 
DWA amplitudes are then 



JfJi 



E 



'_'|(ii-nii) 



2J. + 1 



5 V ~"/ 

X (ji,j2,mi,-m2|/,A^) {Ji,I,Mi,N\Jf,Mf) 



X PsPtA, 



01 



(22) 



In this, {^} = ji,j2,rni,m2,I{N),S,T with 
j2 being the particle and ji the hole in a 
particle-hole specification of the transition. 



Thus, in our DWA evaluations of 
the charge-exchange scattering of interest, 
namely ^^Ge(p, n) and '^^Se{n,p) to 1"*" states 
given by a QRPA model, we have used 

1. SP wave functions: harmonic oscilla- 
tors with oscillator length of 2.09 fm. 
Those are used to specify both the op- 
tical potentials and the reaction ampli- 
tudes. 

2. Optical potentials (to give the distorted 
waves) are formed with the Melbourne 
effective A^A^ interaction at the rele- 
vant incident particle energies. The 
occupancies of the single particle level 
are automatically given in the QRPA 
due to pairing and configuration mix- 
ing and are contained in the OBDME 
in Eq. ([121) . 

3. The same effective interactions are used 
in evaluations of the charge-exchange 
cross sections. 

4. The of Eq. ([12]) are taken as the 
Sjj^j^j=i depending upon which reac- 
tion, (p, n) or {n,p), is described. 



IV. RESULTS 

We have stressed the importance of us- 
ing an appropriate optical model to define 
the distorted wave functions in DWA evalu- 
ations of the charge-exchange scattering. We 
contend that potentials formed using the g- 
folding procedure are such for incident nu- 
cleon energies in the range ~ 40 to ~ 300 
MeV and for nuclei for which, at the mini- 
mum, sensible models of their ground state 
structures can be specified. 

For the mass-76 nuclei we consider specifi- 
cally, very few nucleon scattering results have 
been reported. We have found data for the 
scattering of 22.3 MeV protons from both 
^"^Ge and ^^Se and for 64.5 MeV pro- 
ton scattering from '''^Se (29| (see fig. [2]). For 
nucleon scattering off of these targets, 22.3 
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MeV may be too low an energy to have con- 
fidence that the reaction processes not in- 
cluded in the ^r-folding method, e.g. coupled- 
channel effects, multi-step processes, and the 
like, may have importance. For energies 40 
MeV and higher, such extra processes have 
not been needed to find good replication of 
elastic scattering data with q-folding model 



evaluations 10, 30, 31, 32, 33, yj, 35|, when 



good structure, and appropriate effective A^A^ 
interactions are used, and a proper treatment 
of the Pauli principle is made. 

In this study, we have chosen four incident 
energies at which the transition to the lowest 
1+ state in ^^As is considered. Those energies 
are 45, 65, 120, and 200 MeV. Subsequently, 
for our investigation of transitions to all 23 1"*" 
state excitations, we have used just the two 
largest energies of 120 and 200 MeV. Pro- 
ton elastic scattering data have been taken 
taken for all four energies and from many 
targets. With most cases, (/-folding model 
analyses [lo| gave good reproductions of the 
observations, especially whenever good mod- 
els for the structure of the target were avail- 
able. Besides results given in the review l3] , 
in more recent studies the (/-folding method 
has been used to assess the neutron excess 
distributions 



37, 38, 39 , to 



in nuclei [3f 
compare with non-relativistic and relativistic 
phenomenological Schrodinger equation solu- 



tions l30l. 31 , and to ascertain neutron halo 



or neutron skin characteristics in light mass 



radioactive nuclei [32, 33 



It is a mantra of ^f-folding studies that 
no adjustments to details specified are con- 
sidered post facto. Consequently, (/-folding 
predictions invariably do not yield the qual- 
ity of fit to a data set that may be ob- 
tained by appropriate adjustment of parame- 
ters in current, phenomenological, optical po- 
tentials [40i]. Nonetheless, the (/-folding ap- 
proach does give cross sections that compare 
well with observations; well enough that 
results 



m 



some cases. 



321, 1331, m 



39 



re- 



vealed whether a nucleus had a neutron skin 
or halo. 




FIG. 2: Ratio to Rutherford cross sections for 
64.5 MeV protons scattering from ''^Se. The 
solid line represents the calculation result of this 
work. 



An analysis of cross sections from the elas- 
tic scattering of 65 MeV protons from ^^Se 
is illustrative of the quality of the (/-folding 
potential results. That data [2^, in ratio to 
Rutherford form, are compared in Fig.[2]with 
the (single calculation) result from our g- 
folding optical model potential of the system. 
With the exception of the forward peak, our 
prediction compares favorably with the result 
of the phenomenological optical potential cal- 
culation of Ogino et al |29(|. The agreement 
with data from our non-phenomenological 
approach suffices to give confidence that the 
non-local, complex, optical potential formed 
fully microscopically is a credible, physically 
justified, one. It is important to note that 
there is no addition of any phenomenologi- 
cal elements as used in what r nay be termed 
semi-microscopic methods (25I . 4ll . 
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A. {p, n) reactions to l"*" states in '''^ As 

Differential cross sections evaluated at 
zero-degree scattering and the total reaction 
cross sections from '^^Ge(p, t?,) and ^^Se(n,p) 
leading to the first l"*" state in ''^As are dis- 
played in Fig. [31 The results found at energies 
of 45, 65, 120, and 200 MeV, are connected by 
solid lines (0° cross sections) and by dashed 
lines (reaction cross sections). The results 
of ''^Ge(p, n) to the first excited 1+ state in 
''^As are larger than those of ''^Se(r2,,]?) to the 
same first excited 1"*" state. Over these en- 
ergies, those ratios range from 25 to 65. It 
is intriguing that both the zero-degree dif- 
ferential cross sections, which increase with 
energy, and the reaction cross sections, which 
decrease accordingly, have such similar ratios. 




50 



100 150 
E (MeV) 



200 



FIG. 3: (Color online) Zero-degree differential 
cross sections (filled circles) and total reaction 
cross sections (filled triangles) from ''^Ge(p, n) 
and '^6Se(n,p) leading to ^6As(1+, 0.947 MeV). 



The scale factors are not simply a zero-degree 
phenomenon. This is emphasized by the dif- 
ferential cross-section results for the four en- 
ergies displayed in Fig. H] for a small range of 
momentum transfer gem from 0. The results 



in the top of this figure for the four ener- 
gies as indicated, are those from the charge- 
exchange {p, n) reaction. The other (smaller 
in magnitude) results are differential cross 
sections for the {n,p) reaction. The scale fac- 
tors are an effect of the Pauli principle. In 
32Ge44, the proton Fermi surface lies between 
the 1^3/2 and the O/5/2 levels while the neu- 
tron Fermi surface lies within the Ogg/2 single- 
particle state. For the {p, n) charge-exchange 
reaction, one must move a neutron into a 
proton level. In these nuclei the GT tran- 
sition operator r~cr can make a nucleon into 
a single particle level with the same quan- 
tum numbers or to the spin-orbit partner. 
This is possible for transitions O/5/2 — > O/5/2, 
lPi/2 lpi/2 and Ogg/2 — > OgQ/2- But for the 
inverse reaction, {n,p) on 34Se42, all the possi- 
ble single-particle GT transitions are strongly 
Pauli hindered, if not Pauli blocked. The 
latter cases allow GT transitions since the 
GT operator can only move a proton into 
a corresponding neutron level of the same 
{n£j) orbit or the spin-orbit partner because 
of the smearing of the Fermi surface which is 
mainly induced by pairing correlations. For 
the charge-exchange reaction {n,p) on the 
other hand, the finite range character of the 
transition operator and the knock-out pro- 
cess associated with antisymmetrization al- 
low non-GT type transitions to contribute. 



The energy variation in magnitudes of these 
results for each reaction separately reflects 
the energy dependence of the contributing 
terms to the charge exchange process in the 
effective A^A^ interaction. Tensorial compo- 
nents, which do not contribute strongly to 
the charge exchange processes we consider 
herein, nonetheless vary in importance in the 
overall prescription of the medium modified 
force lOj with energy and with angular mo- 



mentum transfer 34, 35 . Above 100 MeV 



these cross sections are essentially of the same 
magnitude and so we consider the two en- 
ergies, 120 and 200 MeV, in the subsequent 
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Ge(p,n)l 

^200 



10'' 




45 



0.1 



0.2 



FIG. 4: (Color online) Differential cross sections 
from the two charge-exchange reactions leading 
to the first 1"*" state in ^^As for bombarding en- 
ergies of of 45 (solid curves), 65 (dashed curves), 
120 (dot-dashed curves) and 200 (double dot- 
dashed curves) MeV. 



discussions. We note that, for all four ener- 
gies, and for those of 120 and 200 in partic- 
ular, cross sections smoothly decrease away 
from the zero-momentum transfer values over 
the small momentum transfer values consid- 
ered. Furthermore, while the momentum 
transfer values for zero-degree scattering at 
these two energies are small but non-zero, 
evaluations of cross sections setting energies 
to give a zero-momentum transfer gave val- 
ues less than a percent different from those 
of the actual zero-degree calculations. 

The foregoing dealt only with the exci- 
tation of the first 1"^ state in ^^As. We 
now consider the zero-degree cross sections 
for all 23 possible 1"*" states defined by the 
QRPA and their ratios to the correspond- 
ing GT strengths. Shown in table [TTl those 
cross-section values for the ™Ge(p, n) and 
^^Se(ri,jr)) reactions to each of the 23 excited 
l"*" states are listed in columns 2 and 4. The 
ratios of those with the corresponding GT 
strengths are listed in columns 3 and 5. 



TABLE II: differential cross sections at zero 
degrees scattering for the charge-exchange 
reactions '^^Ge{p,n) and '^^Se{n,p) exciting 
^^As(l~'', m). The projectile energy in all cases 
was 200MeV. V are the ratios of each of those 
cross sections with the associated, dimensionless, 
GT strength. The calculated excitation energy 
Ex of the l"*" states is measured from the ground 
state of ^'^Ge. 



m 


Ex, MeV 


^6Ge(p,n)(0°) 


V 


^6Se(n,p)(0°) 


V 


1 


1.16 


5.61 


3.73 


0.09 


3.15 


2 


2.10 


1.74 


5.55 


1 X 10-^ 


25.48 


3 


2.41 


1.01 


5.27 


0.03 


10.96 


4 


2.87 


2.35 


6.19 


0.22 


6.93 


5 


3.14 


1.07 


3.80 


3 X lO"'' 


571.43 


6 


3.41 


3.04 


4.95 


3 X 10-^ 


40.12 


7 


3.95 


12.28 


4.16 


0.26 


3.71 


8 


4.86 


2.31 


3.56 


0.08 


3.54 


9 


5.16 


15.01 


3.96 


0.38 


3.10 


10 


6.39 


0.34 


5.45 


0.18 


6.27 


11 


8.40 


18.43 


3.47 


0.04 


3.40 


12 


9.90 


4.77 


4.09 


0.31 


3.83 


13 


11.25 


3.98 


4.21 


0.24 


4.29 


14 


11.44 


9.40 


4.01 


0.01 


4.77 


15 


12.31 


10.05 


4.05 


0.05 


3.63 


16 


12.60 


51.52 


4.02 


0.01 


3.91 


17 


12.82 


1.49 


3.93 


3 X 10^^ 


8.93 


18 


13.47 


0.08 


3.94 


2 X 10^3 


4.25 


19 


13.63 


0.14 


4.14 


0.01 


4.16 


20 


14.37 


0.27 


3.98 


0.02 


3.92 


21 


15.01 


5 X 10"^ 


4.10 


2 X 10"^ 


3.20 


22 


16.71 


0.02 


4.87 


7 X 10-^ 


6.13 


23 


17.36 


0.10 


4.21 


0.02 


4.58 



State 16 corresponds to excitation of the GT 
resonance that is reflected by the large [p, n) 
cross section at zero-degree scattering being 
51 mb/sr. The corresponding GT strength 
has a dimensionless value of 12.8. The states 
7, 9, 11, 15 and 16 have differential {p,n) 
cross sections at zero-degree scattering larger 
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than 10 mb/sr. The ratio to the GT strength 
for these five states are 4.16, 3.96, 3.47, 4.05 
and 4.02. The ratio for these five leading 
states therefore hes between 3.47 and 4.16; 
and so there is about 20% variation relative 
to the mean value. The six largest values 
for the ^^Se(?T,,p) cross section at zero-degree 
scattering are obtained for the states identi- 
fied as 4, 7, 9, 10, 12 and 13 in the sequence. 
The ratios of these cross sections to their cor- 
responding GT strengths vary between 3.1 
and 6.9. Thus for these six largest {n, p) tran- 
sitions, there is a variation in the ratio of the 
zero-degree charge-exchange cross section to 
the GT strength of about 80%. This is large 
in comparison to the variation in the ratios 
involving the strongest (p, n) reaction cross 
sections. 

In fig. [5] the zero-degree cross sections 
for ''^Ge{p,ny^As (l+,m); m = 1,...,23 
are shown by the filled circles connected by 
dashed lines. The filled squares connected by 
solid lines are the dimensionless GT strengths 
of the operator of eq. ([13]), B{GT^-\ ™Ge 
^™As), for transition to each state of the 
QRPA given in sequence in table [Tll The 
charge-exchange cross section values, con- 
nected by the dashed lines, resulted from 
DWA calculations made using the full NN 
interaction (Melbourne force) as the transi- 
tion operator. The open circles are results 
obtained when only the central part of that 
transition operator was used. Clearly, for 
zero-degree scattering, the two-body spin- 
orbit and tensor contributions do not ef- 
fect the cross sections appreciably. The pro- 
ton incident energy for all (p, n) reactions is 
200MeV. 



Clearly the zero-degree cross sections for all 
(p, n) transitions track similarly to the GT 
strengths of the same states. For the (p, n) , 
and presumably also for the corresponding re- 
actions (^He, t) , the proportionality between 
those charge-exchange reaction cross sections 
in the forward direction and the GT strengths 




QRPA state number 



FIG. 5: (color online) differential cross sections 
at zero-degree scattering from DWA calculations 
of ^^Ge(p, n)^^As(l^, m). Results found with 
and without non-central components in the tran- 
sition operator are depicted by the filled and 
open circles respectively. The filled squares are 
the values of B{GT^-^) for ■^6Ge^'^^As(l+, m) 
for each 1"*" state. 

is fulfihed quite well. As noted above, that 
means a proportionality within about 20% for 
the five strongest transitions. 

The situation is different for the (n,p) 
reactions, and presumably also for the cor- 
responding reactions (t,^He) and (^Li,^Be). 
That is evident both from inspection of the 
results in tablelTTland in fig. O For these tran- 
sitions, most components are Pauli forbidden 
so far as the GT operator is concerned. Fi- 
nite values occur only due to a smearing of 
the Fermi surfaces. But the (n,p) reactions, 
while hindered similarly, also can proceed by 
excitation of other components in the wave 
functions. 



From fig. [6] it is evident that the (n,p) cross 
sections variation over the 23 QRPA possi- 
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5 10 15 20 
QRPA state number 



FIG. 6: (color online) differential cross sections 
at zero-degree scattering for ''^Se(n,p) to the 
23 1^ states in '^^As (filled circles connected by 
dashed lines), the DWBA results found by re- 
stricting the transition operator to just the cen- 
tral terms are depicted by the open circles. The 
associated GT strengths are shown by the filled 
squares connected by solid lines. 

ble states still tracks the values of the associ- 
ated GT strengths. However, note that both 
cross sections and GT strengths are much 
smaller than their counterparts in fig. and 
the omission of non-central force elements in 
DWA calculations makes some greater varia- 
tion than seen with the {p, n) results. Con- 
sequently, the proportionality 'constant' of 
the ratio of forward direction, {n,p) charge- 
exchange cross sections to the GT strengths 
for the six strongest transition is very large. 

As the specifics of the Melbourne force 
change with incident energy, and particularly 
those of the non-central components, we have 
made DWA calculations at other energies. At 
120MeV, the zero-degree cross section val- 
ues and their ratios with the associated GT 
strengths are listed in table lllli As with 
the 200MeV results, the proportionality be- 



tween the charge-exchange cross sections cal- 
culated at zero-degree scattering and the GT 
strength is fulfilled to within about 20% for 
the reaction ^^Ge(p, n). However, the vari- 
ation is much larger for the ratios with the 
cross sections for '^^Se{n,p). In fact the vari- 
ation of the latter ratios is near 90% when 
one considers only the strongest {n,p) tran- 
sitions. 

TABLE III: the zero-degree charge-exchange 
cross sections and ratios (V) of them to the asso- 
ciated GT strength for the ™Ge(p, n) ^^As and 
'^^Se{n,p) ^^As reactions to '''^As(l"'", m). In this 
case, the incident energy was 120MeV. 



m 


^6Ge(p,n)(0°) 


V 


^6Se(n,p)(0°) 


V 


1 


6.51 


4.33 


0.10 


3.58 


2 


2.40 


7.64 


9x10-3 


210.10 


3 


1.47 


7.67 


0.05 


17.26 


4 


2.98 


7.85 


0.28 


8.80 


5 


1.32 


4.70 


6x10-^ 


1146.9 


6 


3.73 


6.07 


2x10-3 


25.56 


7 


14.85 


5.03 


0.31 


4.35 


8 


2.77 


4.26 


0.09 


4.13 


9 


17.47 


4.61 


0.46 


3.73 


10 


0.46 


7.45 


0.24 


8.68 


11 


21.52 


4.06 


0.04 


3.91 


12 


5.65 


4.84 


0.37 


4.56 


13 


4.81 


5.08 


0.24 


4.19 


14 


11.07 


4.73 


0.02 


8.06 


15 


11.83 


4.76 


0.07 


4.50 


16 


61.63 


4.81 


0.02 


6.29 


17 


1.84 


4.86 


2x10-3 


49.43 


18 


0.09 


4.51 


2x10-3 


5.12 


19 


0.17 


5.06 


0.01 


4.68 


20 


0.33 


4.87 


0.02 


4.91 


21 


6x10-3 


4.90 


7x10-5 


10.88 


22 


0.02 


5.68 


9x10-^ 


8.10 


23 


0.11 


4.54 


0.02 


4.65 



In fig. [71 the zero-degree cross sections from 
our DWA evaluations of all '^^Se{n,p) reac- 
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tions to the QRPA 1"*" states in '''^As are 
shown for an incident neutron energy of 
120MeV. As before, those results are depicted 
by the filled and open circles (connected by 
the dashed lines to guide the eye), with the 
filled circle presenting the results when the 
complete Melbourne force is used and the 
open circles giving the results when only the 
central force components are considered. The 
associated GT strengths are depicted by the 
filled squares connected by the solid lines. As 
with the results for 200MeV, the {n,p) cross 
section values vary across the 23 QRPA cases 
very similarly to the GT strength values. But 
the devil is in the differences again and the 
ratio of them is far removed from being con- 
stant over the set. 



^^Ge — >^^As (filled circles connected by a 
solid lines) to the 23 QRPA states in '^^As and 
for the incident proton energy of 200MeV, are 
shown. The ratios for ^^Se(n,p) ^^As cross 
sections are displayed by the open diamonds 
connected by dashed lines (200MeV) and by 
the filled triangles connected by dot-dashed 
lines (120MeV). On this scale the relative 
smoothness of the ratios for all 23 QRPA 
cases of ^^Ge — >™As is apparent. The varia- 
tion though is ~ 20%. But the extreme vari- 
ation over the set for the '^^Se{n,p) ™As ra- 
tios makes it impossible to consider such as 
forming a proportionality constant. Even re- 
stricting consideration to the five strongest 
{n,p) transitions yields a variation of ~ 90%. 




5 10 15 20 
QRPA state number 
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FIG. 7: (color online) zero-degree differential 
cross sections from DWA calculations of 120MeV 
'^^Se{n,p) to 1+ excited states in ™As and the 
associated GT strengths. Details are given in 
the text. 



Finally, in figure [8] the ratios of our calcu- 
lated zero-degree charge-exchange cross sec- 
tions and GT strengths for the transition 
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FIG. 8: (Color online) Ratios of cross sections 
from '^^Ge{p, n) ^^As evaluated at zero degrees 
with the GT strengths (filled circles) for an inci- 
dent energy of 200 MeV and for '^^Se(n,p) '^^As 
displayed for 200 MeV (open diamonds) and for 
120 MeV (filled triangles). 
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V. CONCLUSIONS 

In conclusion, from our analysis of excita- 
tions with the QRPA model structure for the 
three mass-76 nuclei, the proportionality be- 
tween forward differential cross sections and 
GT strengths is fulfilled to within about 20% 
for the ^^Ge(p, n) GT transitions. But this is 
not the case with forward differential charge- 
exchange {n,p) reactions unlike the cases in 
light mass nuclei. In light nuclei, where pro- 
tons and neutrons fill the same or nearly the 
same single-particle levels, the proportional- 
ity of forward charge-exchange reaction cross 
sections and the GT strengths seems valid for 
both (n, p) and {p, n) processes. For them the 
proton and neutron Fermi surfaces are sim- 
ilar. But that is not so for medium heavy 
nuclei such as the mass-76 set we considered 
herein. With them the proton and neutron 
Fermi surfaces are quite different, and the ef- 
fects of Pauli-blocking (hindering for smeared 
surfaces) allows the proportionality to be 
good (within 20%) for (p, n) processes but 
not for the (n,p) transitions. This arises be- 
cause the charge-exchange {n,p) transitions 
are sensitive to radial overlaps of single par- 
ticle wave functions whereas the GT values 



are not, and for the (n,p) processes those in 
which a nucleon stays within a given orbit but 
changes type are not as dominant as with the 
(p, n) ones. 

We anticipate that such will be the case for 
systems that have sizeable neutron excess. It 
will be interesting not only to apply the ap- 
proach we have taken in cases whenever rele- 
vant data are available with which additional 
tests of the quality of the nuclear structure 
model can be made but also in any case when 
application of theories for charge- exchange 
processes initiated by composite projectiles, 
as comparably physically justified, can be 
made. 
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